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ABSTRACT: Coronavirus (CoV) entry is mediated by the viral spike (S) glycoprotein, a class I viral fusion 
protein. During viral and target cell membrane fusion, the heptad repeat (HR) regions of the $2 subunit 
assume a trimer-of-hairpins structure, positioning the fusion peptide in close proximity to the C-terminal 
region of the ectodomain. The formation of this structure appears to drive apposition and subsequent 
fusion of viral and target cell membranes; however, the exact mechanism is unclear. Here, we characterize 
an aromatic amino acid rich region within the ectodomain of the $2 subunit that both partitions into lipid 
membranes and has the capacity to perturb lipid vesicle integrity. Circular dichroism analysis indicated 
that peptides analogous to the aromatic domains of the severe acute respiratory syndrome (SARS)—CoV, 
mouse hepatitis virus (MHV) and the human CoV OC43 S2 subunits, did not have a propensity for a 
defined secondary structure. These peptides strongly partitioned into lipid membranes and induced lipid 
vesicle permeabilization at peptide/lipid ratios of 1:100 in two independent leakage assays. Thus, partitioning 
of the peptides into the lipid interface is sufficient to disorganize membrane integrity. Our study of the $2 
aromatic domain of three CoVs provides supportive evidence for a functional role of this region. We 
propose that, when aligned with the fusion peptide and transmembrane domains during membrane 
apposition, the aromatic domain of the CoV S protein functions to perturb the target cell membrane and 
provides a continuous track of hydrophobic surface, resulting in lipid-membrane fusion and subsequent 


viral nucleocapsid entry. 


The global outbreak of severe acute respiratory syndrome 
(SARS) began in the fall of 2002 in the Guangdong Province 
of China. By July of 2003, 8098 probable cases had been 
reported, resulting in 774 deaths in 29 countries worldwide 
(1). The etiological agent of SARS was quickly identified 
as belonging to the family Coronaviridae (2—7), a group of 
large enveloped RNA viruses exhibiting a broad host range 
and capable of causing respiratory, hepatic, and enteric 
diseases (8). Because of its recent emergence, relative ease 
of transmission, and disease severity, it is necessary to gain 
a better understanding of the pathogenesis of this virus to 
develop antiviral drugs and vaccines for the treatment and 
prevention of SARS. 

All viruses must bind and invade their target cells to 
replicate. Entry of enveloped animal viruses requires fusion 
between the viral membrane and a cellular membrane, either 
the plasma membrane or an internal membrane. Fusion is 


* This work was supported by the National Institutes of Health 
(AI054626, AI054238, RRO18229, and CA08921 to R.F.G. and 
GM60000 to W.C.W.). B.S. is a recipient of a National Research 
Service Award from the NIH (AIO543818). 

* To whom correspondence should be addressed. Telephone: 504- 
988-7076. Fax: 504-988-2739. E-mail: wwimley @tulane.edu. 

* Department of Microbiology and Immunology, Program in Mo- 
lecular Pathogenesis and Immunity, Tulane University Health Sciences 
Center. 

’ Department of Biochemistry, Tulane University Health Sciences 
Center. 

"Louisiana State University Health Sciences Center. 


10.1021/bi0485 15g CCC: $30.25 


mediated by viral fusion proteins, of which two classes have 
been identified for enveloped RNA viruses. Class I fusion 
proteins contain (1) a fusion peptide at or near the amino 
terminus consisting of a stretch of ~20 hydrophobic amino 
acids, (2) a pair of extended a helices, specifically 4,3- 
hydrophobic heptad repeats (HR), and generally, (3) a cluster 
of aromatic amino acids proximal to (4) a hydrophobic 
transmembrane (TM) anchoring domain (9—/2). Gallaher 
and colleagues have shown that several otherwise disparate 
viruses, including orthomyxoviruses, paramyxoviruses, ret- 
roviruses, arenaviruses, and filoviruses, all encode class I 
fusion proteins, varying in length and sequence but highly 
similar in overall structure (/3, /4). The class II fusion 
proteins of alphaviruses and flaviviruses also contain an 
internal fusion peptide but are comprised of three antiparallel 
B-sheet domains (15—17). 


Coronaviruses rely solely on the membrane-bound viral 
spike (S) glycoprotein to bind and infect host cells. The $1 
subunit of the S glycoprotein mediates receptor binding (18, 
19), while the S2 subunit is responsible for driving viral and 
target cell membrane fusion (20). Recent studies indicate that 
the HR regions of the S protein of murine hepatitis virus 
(MHV) and SARS—CoV form a trimer-of-hairpins, a hall- 
mark structure observed for class I viral fusion proteins of 
other RNA viruses. Using peptides corresponding to the two 
HR regions of the $2 subunit, Bosch et al. showed that HR1 
and HR2 of the MHV S2 subunit assemble into a stable 
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antiparallel oligomeric complex (2/). Liu et al. demonstrated 
a similar association between the two HR regions of SARS— 
CoV (22). Most recently, using systematic peptide mapping, 
Tripet et al. showed that the site of interaction between HR1 
and HR2 of the SARS—CoV S2 subunit is between residues 
916—950 of HRI and residues 1151—1185 of HR2 (23). 
They propose that, upon receptor binding, the two HR 
regions of S2 assume a collapsed antiparallel 6-stranded 
a-helical structure (trimer-of-hairpins) that draws the target 
cell and viral membranes together (23). The subsequent 
membrane juxtaposition is believed to result in the destabi- 
lization of both membranes, causing fusion and/or pore 
formation and entry of the viral nucleocapsid into the target 
cell (12, 24-28). These results and additional analyses 
presented herein demonstrate inherent similarities between 
the CoV S protein and fusion proteins of other RNA virus, 
such as HIV (/3, 14, 29), indicating that the CoV S protein 
is a class I viral fusion protein. 

A trimer-of-hairpins structure is a recurring theme for class 
I viral fusion proteins. Trimer-of-hairpins have been identi- 
fied through X-ray crystallographic methods for the filovirus 
Ebola (30—32), oncogenic retroviruses (33, 34), lentiviruses 
HIV and SIV (30, 35, 36), orthomyxovirus influenza (37), 
paramyxovirus SV5 (38), and most recently for MHV (39) 
and SARS—CoV (40). In each case, formation of the trimer- 
of-hairpins brings the N-terminal region of the viral fusion 
protein into close proximity to the C-terminal region of the 
ectodomain, postioning both the fusion peptide and TM 
domains at the same end of the molecule. The apposition of 
both the target cell and viral membranes is believed to result 
in membrane fusion; however, the exact mechanism is 
unclear but likely involves the formation of a fusion pore 
through which the viral nucleocapsid enters the target cell 
(41—46). Several studies (47—50) suggest that the insertion 
of the fusion peptide alone into the target cell membrane 
could potentially cause the membrane distortion necessary 
for fusion, because of its high propensity to partition into 
the membrane interface as measured by the interfacial 
hydrophobicity scale developed by Wimley and White (WW) 
(57). In addition to the fusion peptide, however, a second 
region consisting of aromatic amino acids proximal to the 
TM domain in the fusion proteins of HIV and EboV also 
demonstrate a high tendency to partition into the membrane 
interface according to the WW interfacial hydrophobicity 
scale (5/—54). This region of aromatic amino acids, accord- 
ing to the trimer-of-hairpins model, would also align with 
the fusion peptide and TM domain during apposition of the 
target cell and viral membranes and thus could add to the 
overall hydrophobicity of the environment and contribute to 
the distortion of the lipid membranes necessary for fusion. 
Using synthetic peptides analogous to the aromatic domains 
of HIV-1 gp41 and EboV glycoprotein 2 (GP2), Nieva and 
colleagues have shown that these aromatic regions partition 
into and perturb the integrity of lipid vesicles (10, 55, 56). 
Not surprisingly, proximal to the TM domain of the CoV S 
protein, there is a similar region enriched in aromatic amino 
acids, which is extraordinarily conserved throughout the 
Coronaviridae and lies in an identical location to the 
aromatic domains of HIV (Figures | and 2) and EboV 
(Figure 2). 

Previous studies identifying and characterizing functional 
domains of class I viral fusion proteins of Influenza virus 
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(9, 47, 57), HIV (10, 48, 49, 56, 58), and EboV (50, 55, 59) 
have done so in part using peptides analogous to specific 
regions of these proteins. Because synthetic peptides have 
been shown to be good models to study membrane-associated 
structures and processes functioning in viral fusion (2/—23, 
26, 35, 60-64), we investigated the biochemical and 
functional characteristics of peptides analogous to the 
aromatic domain of SARS—CoV, MHYV, and the human CoV 
OC43. We show that this conserved and highly aromatic 
region of the CoV S2 subunit strongly partitions into the 
membranes of lipid vesicles with a preference for vesicles 
containing anionic lipids, such as phosphatidyl-inositol or 
phosphatidyl-glycerol. Further, we demonstrate that this 
region alone perturbs membrane integrity and causes leakage 
of vesicle contents in two independent leakage assays. Our 
study of the $2 aromatic domain of three CoV is consistent 
with previous findings for HIV and EboV (J0, 55) and 
provides putative evidence for a functional role of this 
aromatic domain in viral fusion and entry. In addition, these 
studies provide additional evidence to support the hypothesis 
set forth by Bosch et al. that the CoV S2 subunit is a class 
I viral fusion protein. 


EXPERIMENTAL PROCEDURES 


Peptide Synthesis. The CoV aromatic (SARSaro, MH Varo, 
and OC43,,0) and SARS—CoV scrambled (SARS,.;) peptides 
were synthesized by solid-phase methodology using a 
semiautomated peptide synthesizer and conventional N-a- 
9-fluorenylmethyloxycarbonyl (Fmoc) chemistry by Gen- 
emed Synthesis, Inc. (San Francisco, CA). Peptides were 
purified by reversed-phase high-performance liquid chro- 
matography, and their purity was confirmed by amino acid 
analysis and electrospray mass spectrometry. Peptide stock 
solutions were prepared in dimethyl sulfoxide (DMSO)! 
(spectroscopy grade), and concentrations were determined 
spectroscopically (SmartSpec 3000, BioRad, Hercules, CA). 


Large Unilamellar Vesicle (LUV) Preparation. LUVs 
consisting of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphati- 
dylcholine (POPC) with 1-palmitoyl-2-oleoyl-sn-glycero-3- 
phosphatidylglycerol (POPG), L-a-phosphatidylinositol from 
bovine liver (PI) (Avanti Polar Lipids, Birmingham, AL), 
and/or cholesterol (CHOL) (Sigma, St. Louis, MI) were 
prepared according to the extrusion method of Nayar et al. 
(65, 66). Briefly, lipids were dried from chloroform solution 
with a nitrogen gas stream and high vacuum overnight. Lipid 
vesicles used in peptide-binding assays and circular dichro- 
ism (CD) experiments were resuspended in 10 mM potassium 
phosphate buffer at pH 7.0 to bring the concentration to 100 
mM total lipid. Samples were subjected to repeated freeze 
and thaw for 15 cycles followed by extrusion through 0.1 
lum polycarbonate membranes in a Lipex Biomembranes 
extruder (Lipex Biomembranes, Vancouver, BC). To prepare 
terbium(III) chloride hexahydrate (Tb**) LUVs, lipids were 
resuspended to 100 mM concentration in 50 mM Tb?*, 100 


' Abbreviations: ANTS, disodium 8-aminonaphthalene-1,3,6-trisul- 
fonic acid; DMSO, dimethyl] sulfoxide; DPA, 2,6-pyridine dicarboxylic 
acid; DPX, p-xylene-bispyridinium bromide; PI, L-a-phosphatidyli- 
nositol from bovine liver; POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3- 
phosphatidylcholine; POPG, 1-palmitoyl-2-oleoyl-sn-glycero-3-phos- 
phatidylglycerol; Tb**+, terbium(III) chloride hexahydrate; TES, 
N-tris(hydroxymethyl)methyl-2-aminoethane sulfonic acid. 
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FiGurE 1: Schematic comparison of the HIV-1 TM with the SARS—CoV S protein. (Left) Updated model of HIV-1 TM from Gallaher et 
al. and Garry et al. (/3, 17). (Right) Hypothetical model of SARS—CoV S (S2 subunit) showing motifs shared with HIV-1 TM. The amino 
terminus of the S2 portion is not depicted. The two HR a-helical regions, N helix (HR1, orange) and C helix (HR2, yellow), are depicted 
at amino acid residues 902—1011 and 1131—1185, respectively. This is consistent with the HR predictions of both Tripet et al. (23) and 
Liu et al. (22). An interhelical domain of approximately 130 amino acids is depicted between the N and C helices. This region is extremely 
similar to the interhelical region for retrovirus TM proteins and EboV GP2 and has therefore been modeled as a similar disulfide-stabilized 
apex. Just prior to the TM anchor (indigo) of S2, there is a region enriched in aromatic amino acids. This region, termed the aromatic 
domain (green), is highly conserved throughout the Coronaviridae and lies in an identical location to the aromatic domains of HIV and 
EboV. The S1 subunit, which includes the receptor-binding domain, is depicted schematically as a large ellipse, corresponding to the 
characteristic large globular headgroups seen in electron micrographs of CoV. 


mM sodium citrate, and 10 mM J-tris(hydroxymethyl)- 
methyl-2-aminoethane sulfonic acid (TES) at pH 7.2. Gel 
filtration on Sephadex G-200 was used to remove unencap- 
sulated terbium in a buffer of 10 mM TES and 325 mM 
NaCl (67). Disodium 8-aminonaphthalene-1,3,6-trisulfonic 
acid (ANTS)/p-xylene-bispyridinium bromide (DPX) LUVs 
were prepared in an analogous manner using a resuspension 
buffer of 25 mM ANTS and 5 mM DPX (Molecular Probes, 
Eugene, OR) in 10 mM potassium phosphate at pH 7.0 (68). 
LUVs were eluted from a Sephadex G-200 gel column using 


10 mM potassium phosphate at pH 7.0. Final lipid concen- 
trations were determined by phosphate analysis (69, 70). 
Peptide-Binding Assay. Partitioning of peptides into lipid 
bilayer was monitored by the fluorescence enhancement of 
tryptophan (7/). Fluorescence was recorded at excitation and 
emission wavelengths of 280 and 340 nm, respectively, and 
8 nm bandwidths using an SML Aminco 8100 spectrofluo- 
rometer (Rochester, NY). Quartz cuvettes were used with 
excitation and emission path lengths of 4 and 10 mm. The 
measurements reported here were carried out in 10 mM 
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FiGuRE 2: Hydropathy plots corresponding to sequences of SARS— 
CoV S82, HIV-1 gp41, and EboV GP2. WW interfacial hydropho- 
bicity plots (—) (SJ) and WW octanol hydrophobicity plots (- - -) 
(73) were generated for individual residues of (A) the SARS—CoV 
strain Urbani S2 subunit (amino acids 775— 1225), (B) HIV-1 strain 
HXB2 gp41 (amino acids 502—710), and (C) EboV strain Zaire 
GP2 (amino acids 520—676). Mean values for a window of 13 
residues are plotted. Areas shaded in gray represent aromatic 
domains of the represented class I viral fusion proteins. A = 
aromatic domain, and TM = transmembrane domain. 
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potassium phosphate at pH 7.0. In other buffer systems, the 
coefficients are very similar, indicating that partitioning is 
not sensitive to the exact buffer conditions or ionic strength. 
Peptides were added from DMSO stock solutions to 250 uL 
of buffer and mixed by inversion. Likewise, LUVs were 
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titrated to a final lipid concentration of 1 mM and mixed by 
inversion. Intensity values (J) were adjusted for lipid scat- 
tering and normalized to the peptide in buffer (Jo). Partition- 
ing coefficients were obtained using eq | where K, is a mole 


T/T, = 1+ ((K, x [L/W] + (K, x [ED))) x 
(Cadi 1) 


fraction partition coefficient that represents the amount of 
peptide in bilayers as a fraction of the total peptide present 
in the system, Imax is a variable value for the fluorescence 
enhancement at complete partitioning determined by fitting 
the equation to the experimental data, [ZL] is the concentration 
of the lipid, and [W] is the concentration of water (55.3 M). 
Tb?*/2,6-Pyridine Dicarboxylic Acid (DPA) Microwell 
Plate Assay. For the microwell plate assay, a 200 wL aliquot 
of vesicle solution containing 500 uM Tb?* LUVs in 10 mM 
TES, 50 uM DPA, and 325 mM NaC] at pH 7.2 was pipetted 
into each well of a plastic 8 x 12 format plate (67). Peptides 
in DMSO were added to each well at peptide/lipid (P/L) 
molar ratios of 1:500, 1:250, 1:100, and 1:50; well contents 
were thoroughly mixed, and plates were allowed to incubate 
at room temperature for 2 h. In addition to peptide-treated 
wells, DMSO-treated and Triton X-100-treated (Sigma) wells 
served as negative and positive controls, respectively. After 
incubation for 2 h, Tb**/DPA fluorescence was visualized 
under horizontally mounted short-wave (254 nm) ultraviolet 
(UV) light sources in a darkroom (67). Plates were photo- 
graphed, and images were recorded on a Nikon Coolpix 995 
using a 4 s exposure time with 100 speed, 2.6 aperture, and 
a 540 nm band-pass optical filter between the sample and 
lens. For each experimental plate, the Tb>*/DPA fluorescence 
for the peptide-treated wells was compared to that of wells 
containing the same amount of untreated vesicles and to wells 
containing vesicles that had been lysed with the detergent 
Triton X-100. Color adjustment and contrasting were nor- 
malized to negative controls using Adobe Photoshop. 
ANTS/DPX Leakage Assay. ANTS and DPX were en- 
trapped in LUVs at concentrations of 25 and 5 mM, 
respectively. Peptides were added to 500 uM lipid solution 
in 10 mM potassium phosphate at pH 7.2 and P/L molar 
ratios of 1:500, 1:250, 1:100, and 1:50. Although most 
ANTS/DPX leakage occurred in the first several hours after 
peptide addition, samples were shaken overnight at room 
temperature to allow the leakage to go to completion. 
Fluorescence was recorded at excitation and emission 
wavelengths of 350 and 510 nm, respectively, and 8 nm 
bandwidths using an SML Aminco 8100 spectrofluorometer. 
Quartz cuvettes were used with excitation and emission path 
lengths of 4 and 10 mm, respectively. The percent change 
in fluorescence was defined by eq 2 where Fynax is obtained 


% fluorescence = ((F — Fo)/(Finax ~ Fo) X 100% (2) 


ax 


by adding 25 wL of 5% Triton X-100 and Fo is equivalent 
to DMSO controls. Percent leakage was calculated from the 
fraction released (four), where four Was determined using curve 
fitting by numerically solving eq 3 where fou is the fractional 


% fluorescence = 100 x (fy + C1 — four) X 
1({DPX];,)) (3) 
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leakage from the vesicles and [([DPX]j,) is the fluorescence 
intensity of the fraction of ANTS that remains entrapped. 
The entrapped ANTS is quenched by the DPX that remains 
entrapped, [DPX]in. For each point, [DPX]jn is calculated by 
(1 — four) x [DPX]o, where [DPX ]o is the initial concentration 
of DPX. 


CD Spectroscopy. CD spectra were recorded on a Jasco 
J-810 spectropolarimeter (Jasco Inc., Easton, MD), using a 
1 mm path length, 1 nm bandwidth, 16 s response time, and 
a scan speed of 10 nm/min. All CD runs were performed at 
room temperature with the peptide dissolved in 10 mM 
potassium phosphate buffer at pH 7.0. LUVs were added at 
a lipid concentration of | mM from a stock in 10 mM 
potassium phosphate buffer at pH 7.0 and allowed to incubate 
for ~30 min. Three successive scans between 190 and 250 
nm were collected, and the CD data are expressed as the 
mean residue ellipticity, derived from the formula © = (deg 
x cm’)/dmol. 


Proteomics Computational Methods. Methods to derive 
general models of surface glycoproteins have been described 
previously (/3). Domains with significant propensity to form 
TM helices were identified with TMpred (ExPASy, Swiss 
Institute of Bioinformatics) and Membrane Protein 
eXplorer (MPeX, Stephen White laboratory, http:// 
blanco. biomol.uci.edu/mpex). TMpred is based on a statisti- 
cal analysis of TMbase, a database of naturally occurring 
TM glycoproteins (72), while MPeX detection of membrane- 
spanning sequences is based on experimentally determined 
hydrophobicity scales (5/7, 54). Sequences with a propensity 
to partition into the lipid bilayer were also identified with 
MPeX using interfacial settings. 


RESULTS 


Interfacial Hydrophobicity Analysis and Peptide Se- 
quences. WW hydrophobicity scales were used to identify 
regions of the S2 subunit of the CoV S protein with high 
propensity to partition into lipid membranes. These scales 
are based on the free energies of transfer AG (kcal/mol) of 
amino acid sequences from water into bilayer interfaces and 
into n-octanol, taking into consideration the contribution from 
the peptide bond (5/). The octanol scale is used as a predictor 
of TM domains (73), while the interfacial scale predicts 
interfacial binding (5/). Because of the salient similarities 
between the S2 subunit and the class I fusion proteins of 
other RNA viruses such as HIV-1 (Figure 1), we compared 
the hydrophobicity plots of the SARS—CoV S2 subunit to 
that of the HIV-1 gp41 and EboV GP2 (Figure 2). 


When the average interfacial hydrophobicity was plotted 
for the fusion proteins of these three viruses, similar regions 
with high propensity for membrane partitioning were detected 
(Figure 2). At the N-terminal region of all three fusion 
proteins, a region of high interfacial hydrophobicity was 
detected. For HIV-1 and EboV, this region corresponds to 
the viral fusion peptide (parts B and C of Figure 2). Although 
no putative fusion peptide has been determined for the 
SARS—CoV fusion protein, a stretch of 19 hydrophobic 
amino acids comprised of residues 770 MWKTPTLKYFG- 
GFNFSQIL 788 with an average interfacial hydrophobicity 
score of 3.07 kcal/mol was detected at the N-terminal region 
of S2 (Figure 2A). The location of this region is 9 amino 
acids downstream of the minimum furin cleavage site 
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Table 1: Amino Acid Sequences and WW Hydrophobicity Scores 
of CoV Aromatic Peptides 


net position WWIH? 
peptide aminoacid sequence’ charge withinS2 (kcal/mol) 
SARSaro KYEQYIKWPWYVW +1 1187-1199 3.58 


MHVar0 TYEMYVKWPWYVW 0 
OC43a0 TYEYYVKWPWYVW 0 1289-1301 5.57 
SARS. YEWKWIYWYPVKQ +1 1187-1199 3.58 


“The SARS aromatic (SARSarw), MHV aromatic (MHV,,.), and 
OC43 aromatic (OC43.:0) peptides were synthesized based on their 
amino acid sequence determined from GenBank accession number 
AY278741 (SARS—CoV strain Urbani), AY497331 (MHV strain A59), 
and NP_937950 (human CoV OC43), respectively. Amino acid 
differences among the three CoV aromatic peptides are shown in bold 
and underlined text. ’ Hydrophobicity scores were determined according 
to the Wimley and White interfacial hydrophobicity (WWIH) scale 
using a window of 13 residues. ©“ The SARS,;. sequence was arbitrarily 
scrambled to generate the peptide SARS cer. 


1264—1276 4.86 


(residues 758—761) and almost coincident with that of the 
HIV-1 fusion peptide; therefore, it should be considered as 
a possible fusion domain of the SARS—CoV S protein. 

A second region of high hydrophobicity was detected at 
the C-terminal end of the fusion proteins, correlating to the 
putative TM domain of the SARS—CoV S2 (residues 1190— 
1225 of Figure 2A) and the experimentally determined TM 
anchors of HIV-1 gp41 and EboV GP2 (residues 665—700 
of Figure 2B and residues 644—672 of Figure 2C, respec- 
tively). Nieva and colleagues have shown that, for HIV-1 
and EboV, this large region of high interfacial hydrophobicity 
is segmented into two independent domains: one aromatic 
amino acid rich segment at the C-terminal portion of this 
domain and an adjacent segment comprising the TM anchor 
of the fusion protein (0, 55, 56). The separation of this large 
region into two distinct domains is consistent with the 
observation that the C-terminal hydrophobic domains of these 
viral fusion proteins are considerably longer (35—40 amino 
acids) than the ~20 amino acids required for a single & helix 
to span a membrane. Likewise, WW interfacial and octanol 
hydrophobicity plots of the C-terminal end of the SARS— 
CoV S2 subunit shows a pattern similar to that of the HIV-1 
gp41 and EboV GP2 (part A versus parts B and C of Figure 
2). A 12 amino acid region rich in aromatic amino acids, 
with an average WW interfacial hydrophobicity score of 3.58 
kcal/mol, is present. This region is proximal to a predicted 
TM domain (Figure 2A), which was identified using the WW 
octanol hydrophobicity scale that allows for the identification 
of TM helices of membrane proteins with very high accuracy 
(73). On the basis of these observations, it is unlikely that 
the aromatic domain is part of the TM anchor as previously 
predicted by Rota et al. (6), rather, as described for the 
aromatic domains of HIV-1 and EboV (J0, 55, 56, 58), this 
region is most likely an independent domain proximal to the 
TM anchor of the S2 subunit. 

Sequence analysis of the S2 subunit of MHV and the 
human CoV OC43 showed coinciding interfacial hydropho- 
bicity plots to that of the SARS—CoV S2 subunit (data not 
shown). In addition, the presence of a highly conserved 
aromatic domain, differing in only 3 amino acids to the 
SARS aromatic domain, was identified (Table 1). Interfacial 
hydrophobicity scores of 3.58, 4.86, and 5.57 kcal/mol were 
predicted for the aromatic domains of SARS—CoV, MHV, 
and OC43, respectively. On the basis of the these analyses, 
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FiGuRE 3: SARSzro peptide partitions into membranes of LUVs. 
Change in tryptophan fluorescence of SARSaro peptide as a function 
of increasing concentrations of LUVs composed of (ll) POPC, (@) 
POPC/PI (9:1), (4) POPC/POPG (9:1), or (0) POPC/PI/CHOL (6.5: 
1:2.5). LUVs were titrated at concentrations of 100, 250, 500, 750, 
and 1000 “M lipid with 2.5 uM peptide. Tryptophan fluorescence 
values at each lipid titration (F) were normalized to tryptophan 
fluorescence values in 10 mM potassium phosphate buffer alone 
(Fo). 


peptides of 13 amino acids in length were synthesized and 
used throughout this study to determine the functional 
importance of this region within the CoV S2 subunit (Table 
1). 

CoV Aromatic Domains Interact with Lipid Membranes. 
On the basis of the interfacial hydrophobicity scores of the 
CoV aromatic peptides, we believed that they could partition 
into membranes. We first assessed the ability of the CoV 
aromatic peptides to interact with membranes of LUVs 
composed of different lipid compositions. LUVs composed 
of POPC with PI, POPG, and/or CHOL were used as targets 
in partitioning experiments with the CoV aromatic peptides. 
The degree to which a peptide partitions into a lipid vesicle 
can be determined fluorometrically by observing the change 
in tryptophan fluorescence (F) as a function of the increasing 
lipid concentration. The fluorescence of tryptophan increases 
in a low-polarity environment such as the lipid membrane 
interface. Figure 3 shows the normalized tryptophan fluo- 
rescence (F/F) for the SARS,,. peptide as a function of the 
increasing lipid concentration of different LUVs. SARSaro 
fluorescence exponentially increased as a direct function of 
increasing lipid concentrations. 

The partition coefficients derived from tryptophan fluo- 
rescence curves for each CoV peptide are summarized in 
Table 2. All four CoV peptides partitioned into lipid 
membranes, consistent with their hydrophobicity. The degree 
of partitioning for all four peptides was similar, and the 
presence of anionic lipids (i.e., PI or POPG) enhanced 
peptide partitioning. The addition of CHOL, a natural 
component of mammalian membranes, decreased peptide 
partitioning, most notably with POPC/CHOL LUVs. This 
effect, however, could be partially reversed by the addition 
of small amounts of anionic lipids to the vesicles (Table 2), 
suggesting an intrinsic role for anionic lipids as part of the 
membrane composition. 

All of the structural and functional experiments described 
below were carried out at lipid concentrations between 0.5 
and 1 mM lipid. The partition coefficient measurements 
(Table 2) were carried out at slightly lower peptide concen- 
trations than most of the other experiments. Thus, the 
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measured partition coefficients represent a minimum possible 
value because peptide self-association at higher concentra- 
tions can increase partition coefficients. Partitioning cannot 
be lower than that measured at an equal or lower peptide 
concentration. Because of this, the minimum fraction of 
peptide bound can be calculated from the data in Figure 3 
and Table 2 by fraction bound = (K, x [L])/CK,y x [LZ] + 
[W]), where [LZ] is the lipid concentration and [W] is the 
concentration of pure water (55.3 M). For most of the 
experiments in this work, the minimum fractional binding 
is between 30 and 90%. Fractional binding to CHOL- 
containing POPC was not strong enough to be accurately 
measured and therefore was less than 20%. For clarity, we 
report total system composition when giving P/L ratios. The 
reader should recognize that the actual P/L ratio in the 
membranes will always be lower, especially for CHOL 
containing bilayers. 

Tb?*/DPA Microwell Assay. To test the potential of the 
CoV aromatic peptides to perturb membrane integrity, a high- 
throughput leakage assay was used. The Tb**/DPA microwell 
assay is a stringent, sensitive visual screening assay devel- 
oped to rapidly identify peptides capable of permeabilizing 
lipid membranes (67). The detectability is based on the strong 
fluorescence emission of the lanthanide metal Tb** when it 
interacts with the aromatic chelator DPA. In this binary 
system, Tb** is entrapped in lipid vesicles, while DPA is 
present on the outside. Membrane permeabilization is 
required for the entrapped Tb** to interact with DPA, 
resulting in fluorescence emission. In the experimental assay, 
CoV aromatic peptides were incubated with Tb** LUVs (500 
UM lipid) at P/L molar ratios of 1:100 and 1:50. After 
incubation for 2 h at room temperature, the extent of Tb** 
leakage from lipid vesicles was visually determined by the 
detection of a bright green fluorescence upon irradiation with 
UV light. An example plate is shown in Figure 4 in which 
the SARS.1o (rows 1 and 2) and SARS. (rows 3 and 4) 
peptides were tested for their potential to permeabilize LUVs 
composed of POPC, POPC/PI (9:1), or POPC/POPG (9:1). 
The SARSao peptide at P/L ratios of 1:100 and 1:50 
permeabilized all three LUVs tested, with the greatest degree 
of fluorescence detected in wells with POPC or POPC/PI 
(9:1) LUVs. The extent of leakage induced by SARS,,,. was 
less than the observed leakage in the detergent-solubilized 
wells (row 6 of Figure 4). However, we have found that the 
visual Tb**/DPA assay is often more stringent than typical 
fluorometric assays of membrane permeabilization. Thus, any 
visible response, such as that observed in Figure 4 is taken 
to be a positive result. In contrast, the SARS... peptide did 
not induce leakage of any of the three LUVs tested, as 
detectable by this assay. Comparable results were achieved 
with the MHVaro and OC43.;5 peptides at P/L ratios of 1:100 
and 1:50, with OC43,,. exhibiting slightly lower levels of 
leakage (data not shown). 

ANTS/DPX Leakage Assay. We used the ANTS/DPX 
leakage assay, as described in the Experimental Procedures, 
as a means of quantitating the membrane permeabilization 
capacity of the CoV aromatic peptides. The ability of the 
CoV aromatic peptides (Figure 5) to release the encapsulated 
fluorescent probe ANTS from LUVs was examined at P/L 
ratios of 1:500, 1:250, 1:100, and 1:50. Figure 5 illustrates 
the extent of leakage observed for all four peptides and the 
effect of different lipid compositions on membrane perme- 
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Table 2: Partition Coefficients of CoV Aromatic Peptides 


partition coefficients (x 10%)" 


LUV compositions?’ SARSaro MHVaro OC43.ar0 SARS er 
POPC 49 + 13 101+ 13 52+ 14 47+ 11 
POPC/PI (9:1) 112+ 18 289 + 27 104+ 17 82+ 18 
POPC/POPG (9:1) 142 + 13 367 + 48 21+9 205 + 15 
POPC/CHOL (7.5:2.5) <20 <20 <20 <20 
POPC/PI/CHOL (6.5:1:2.5) 80 + 273 + 36 75 + 44 89 +9 
POPC/POPG/CHOL (6.5:1:2.5) 56 + 22 115+9 26 + 26 49+ 11 


“ Partition coefficients based on spectroscopic measurements of the change in tryptophan fluorescence as a function of increasing lipid titrations. 
Tryptophan intensity values (J) were fitted to eq | and the resulting partition coefficients are presented as the means + standards errors (n = 5 per 
group). ’ Lipid was titrated (100 uM — 1000 uM) with 2.5 uM peptide in a volume of 250 wL. 


—_— 


A B C 
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FiGuRE 4: SARSaro peptide induces leakage of LUVs. Each well 
contained 250 wL of 50 uM DPA and 500 uM Tb**-encapsulated 
LUVs composed of (A) POPC, (B) POPC/PI (9:1), or (C) POPC/ 
POPG (9:1). Wells were treated with SARS, peptide at P/L molar 
ratios of 1:100 or 1:50 (1 and 2), SARS, peptide at P/L molar 
ratios of 1:100 or 1:50 (3 and 4), 20 uL of DMSO (5), or 20 uL of 
Triton X-100 (6). Plates were incubated for 2 h at room temperature, 
and membrane permeabilization was determined by visual detection 
of Tb?+/DPA fluorescence. 


abilization. As with the Tb**/DPA microwell assay, the 
SARSaro peptide induced leakage of ANTS from LUVs to a 
greater degree than its scrambled counterpart, SARS... On 
average, the percent leakage detected at all P/L ratios was 
approximately 2 times greater for the SARS,,. peptide as 
compared to the SARS,., peptide (part A versus part B of 
Figure 5). The degree of leakage induced by SARSaro varied 
based on the lipid composition of the LUVs tested. The 
percent leakage detected from LUVs composed of either 
POPC/PI or POPC/POPG was 34 and 36%, respectively, as 
compared to 29% leakage observed in POPC LUVs at P/L 


07——A SS 0 
0.000 0.005 0.010 0.015 0.020 0,005 0.010 0.015 0.020 
[Peptide:Lipid] (uM) [Peptide:Lipid] (uM) 


FiGure 5: Extent of leakage from ANTS/DPX LUVs induced by 
the CoV aromatic peptides. (A) SARSaro, (B) SARSger, (C) MH Varo, 
or (D) OC43,,. peptides were added to 500 “M LUVs composed 
of (Ml) POPC, (@) POPC/PI (9:1), (&) POPC/POPG (9:1), (G) 
POPC/CHOL (7.5:2.5), (O) POPC/PI/CHOL (6.5:1:2.5), or (A) 
POPC/POPG/CHOL (6.5:1:2.5) at P/L molar ratios of 1:500, 1:250, 
1:100, or 1:50. Samples were incubated at room temperature for 
24 h before measuring the extent of leakage fluorometrically. 
Percent leakage was calculated using eq 3 (n = 3 per group). 
ratios of 1:50 (Figure 5A). In addition, the presence of CHOL 
reduced the amount of leakage induced by the SARSaro 
peptide. As compared to the SARS,,. peptide, MHVaro 
induced equal or greater leakage of ANTS from LUVs, and 
the extent of leakage also appeared to be dependent upon 
the lipid composition of the vesicles. For example, at P/L 
ratios of 1:50, the percent leakage detected in LUVs 
composed of either POPC/PI or POPC/POPG was 50 and 
60%, respectively, as compared to 20% leakage observed in 
LUVs composed solely of POPC (Figure 5C). Leakage in 
the presence of CHOL was reduced for all peptides but 
generally not to the extent predicted by the much lower 
partitioning (Table 2), indicating that the membrane-perme- 
abilizing activity per bound peptide was equal or higher in 
the presence of cholesterol. 

The activity of OC43,;0 in the ANTS/DPX leakage assay 
was considerably less in this assay than for the other two 
CoV aromatic peptides (Figure 5D). At P/L ratios of 1:50, 
OC43.a10 induced minimal ANTS leakage in vesicles without 
cholesterol (~15—20%) and reduced leakage in vesicles 
containing cholesterol (~5—15%). This is in contrast to the 
observation, discussed above, that OC43,,. caused leakage 
when assayed by the Tb**/DPA system. Although there is 
only a single amino acid difference between MHV,,,. and 
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FiGuRE 6: CD spectra (mean residue ellipticity ©) of the CoV aromatic peptides: (A) SARSaro, (B) MHVaro, (C) OC43 aro, and (D) SARS cr 
in 10 mM potassium phosphate buffer at pH 7.0 alone (MM) or with 1 mM LUVs (Q) composed of POPC/PI (9:1) at room temperature (P/L 


ratio of 1:20). 


OC43 10, the pore formed by OC43,,. peptide is apparently 
not the same. Nonetheless, by demonstrating OC43ar0- 
induced leakage in at least one system, our data support the 
conclusion that OC43.0, like the other aromatic domain 
peptides, binds and destabilizes lipid bilayer membranes. 
CD Analysis of the CoV Aromatic Peptides. To examine 
the potential for the formation of secondary structure upon 
interaction with lipid membranes, the CoV aromatic peptides 
were examined by CD spectroscopy. Parts A—C of Figure 
6 show representative far-UV CD spectra of the CoV 
aromatic peptides in buffer and with LUVs, and Figure 6D 
shows the results for the SARS. peptide. In 10 mM 
potassium phosphate buffer at pH 7.0, peptides showed a 
random coil spectrum with single minima at 200 nm. We 
next analyzed the potential of the CoV aromatic peptides to 
adopt a secondary structure in the presence of lipids by 
adding LUVs composed of POPC/PI at lipid concentrations 
of 1 mM. The SARS,,. peptide showed very little change in 
the CD spectrum, despite the fact that at least 60% of the 
peptide is bound to the vesicles under these conditions (Table 
2). Similarly, for the MHVa. and OC43.;. peptides, there 
was only a small change in the observed CD spectra as 
compared to the buffer alone (parts B and C of Figure 6). 
Although not indicative of a predominant, well-defined 
secondary structure, which would have distinct minima at 
208 and 222 nm if @ helical or at ~218 nm if / sheet, it 
appears that the peptides may be assuming a partially ordered 
structure. Nonetheless, the lack of a distinct maximum in 
the CD spectra at ~200 nm demonstrates that these peptides 
are predominantly random coil in membranes. The same is 
true for the SARS,.; peptide, in Figure 6D. The small amount 
of secondary structure apparent in the buffer spectrum does 


not increase when vesicles are added, despite the fact that 
the SARS, peptide binds well to these vesicles. 


DISCUSSION 


The recent emergence of the SARS—CoV, which induces 
a severe respiratory illness with a mortality rate approaching 
10%, has inspired increased interest in all aspects of CoV 
pathobiology. A detailed understanding of mechanisms by 
which CoV S proteins mediate fusion between the virion 
envelope and cellular membrane could contribute to the 
design of antiviral drugs and vaccines. Prior studies suggested 
that CoV S proteins might represent a unique class of viral 
fusion proteins. More recent analyses, however, indicate that 
the CoV S proteins share structural features with class I 
fusion proteins (/2), including a pair of extended a helices 
and a cysteine-linked “hinge” region between the helices (2/, 
23). In addition to these studies, our results demonstrate 
structural and functional similarities between the aromatic 
domain of the SARS—CoV S protein and that of other class 
I viral fusion proteins. Similar to Influenza virus, HIV, and 
EboV, these conserved structural motifs within the SARS— 
CoV S protein presumably play an important role in SARS 
CoV infection, entry, and pathogenesis. 


The 13-residue CoV aromatic peptides studied were 
selected using WW interfacial hydrophobicity plots presented 
in Figure 2. As compared to HIV and EboV, the SARS— 
CoV S2 subunit exhibits common domains with high 
propensities to partition into lipid membranes, of which the 
aromatic domain scores highest for all three viruses. When 
determining the amino acid sequence of the SARS—CoV 
TM domain, Rota et al. described the aromatic domain as 
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part of the TM anchor (6). In Figure 1, we depict the aromatic 
domain as separate from the TM anchor based on WW 
interfacial and octanol hydrophobicity plots of the C-terminal 
end of the SARS—CoV S2 subunit (Figure 2A). The separate 
and distinct peak (6.63 kcal/mol) generated with the WW 
octanol hydrophobicity scale, with a core region correspond- 
ing to amino acids 1200—1217, strongly suggests a TM helix 
adjacent to the aromatic domain. Moreover, 18 amino acids 
are sufficient to constitute a TM anchor as shown in Figure 
1 (54, 74), while the 30 amino acid length of the putative 
TM anchor and the aromatic domain together is considerably 
longer than known single-span TM helices. Therefore, on 
the basis of the relative hydrophobicities (interfacial versus 
octanol) of these two adjacent domains and the total length 
of the C-terminal hydrophobic domain, the SARS—CoV 
aromatic domain most likely constitutes a separate domain 
similar to that of HIV and EboV (J0, 55, 56). 

The aromatic domains studied represent a common motif 
present in some but not all class I viral fusion proteins, which 
is adjacent to or contiguous with the TM anchor domain. 
Peptides corresponding to the aromatic domains of HIV gp41 
and EboV GP? can interact with bilayer membranes and have 
the propensity to disrupt them (J0, 75). Mutational studies 
have demonstrated that viral infectivity is sensitive to the 
sequence and composition of this aromatic domain (75, 76), 
suggesting that the aromatic domain is an important part of 
the virion/cell fusion machinery. In the work presented here, 
we provide data to support the hypothesis that, similar to 
HIV and EboV, the aromatic sequence of the CoV S protein 
also interacts intimately with lipid membranes, induces 
membrane disruption, and may likewise participate in viral 
infectivity. Peptide-binding assays demonstrated strong 
partitioning into lipid membranes composed of the zwitte- 
rionic lipid POPC. Despite their modest net charges of either 
0 (MHV and OC43) or +1 (SARS), the presence of anionic 
lipids (i.e., PI or POPG) increased binding significantly. The 
presence of CHOL, as observed for other viral peptides (55), 
reduced but did not eliminate peptide binding. 

The SARS, MHV, and OC43 aromatic peptides also 
induced measurable lipid vesicle permeabilization in the 
Tb**/DPA microwell assay and varying levels of leakage in 
the ANTS/DPX assay, whereas a scrambled SARS aromatic 
peptide (SARS,.,) induced reduced leakage in both assays. 
It is not surprising that the SARS, aromatic peptide 
interacted with membranes and induced some degree of 
leakage, given its high interfacial hydrophobicity. In fact, 
peptides rich in aromatic amino acids commonly partition 
into membranes as a result of their intrinsic hydrophobicity 
(52). Nonetheless, the results obtained with the SARS,.; 
aromatic peptide suggest that the specific sequence may 
influence the functionality of these peptides in membranes. 
These structural/biophysical studies are consistent with an 
important role for this aromatic domain sequence in CoV- 
mediated fusion and further verify that the CoV S2 protein 
is a class I fusion protein. 

Lipid vesicles containing cholesterol, a major component 
of mammalian membranes in vivo, were used throughout 
the studies described above. Although the addition of 
cholesterol reduced peptide-mediated leakage, as compared 
to cholesterol-free vesicles, measurable leakage in the 
presence of cholesterol was still observed, indicating that 
the membrane-permeabilizing activity per peptide bound is 
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very high in the presence of cholesterol. Several studies have 
examined the effect of cholesterol incorporation in lipid 
vesicles on peptide-mediated leakage of vesicle contents. 
Although some studies report that inclusion of cholesterol 
in membranes reduces leakage (55, 77), other studies report 
increased activity (7S). The latter is observed when choles- 
terol is included along with other lipid molecules, especially 
sphingomyelin. It has been proposed that the addition of these 
two molecules to artificial lipid membranes may mimic the 
formation of segregated membrane phases or lipid “rafts” 
(79-81). Lipid rafts are highly ordered membrane structures 
believed to play a role in receptor-mediated cell signaling 
and viral entry. Although cholesterol is an important 
component of mammalian membranes, it remains to be 
determined how to adequately mimic the complex lipid and 
protein composition of biological membranes in vitro. Studies 
are currently being conducted to investigate the effect of 
sphingomyelin and cholesterol incorporation to the lipid 
vesicle model system described herein. 

We also used CD spectroscopy to show that the coronavi- 
rus aromatic peptides that we studied were not highly 
structured in the presence of membranes, even under 
conditions where the peptides induce significant membrane 
destabilization and leakage. Thus, o-helical secondary struc- 
ture is neither necessary for binding nor necessary for 
membrane permeabilization by these aromatic domains. 
Studies of longer viral aromatic sequences, including peptides 
that contain portions of the amphipathic C helix (HR2), report 
a-helical secondary structure when the peptides are membrane- 
bound (J0, 56, 58). However, secondary structure is a 
strongly driven thermodynamic consequence of membrane 
partitioning for long hydrophobic peptides (54); therefore, 
the biological relevance of observed secondary structure for 
peptides containing the aromatic domain/C helix is not clear. 
It is also not clear what secondary structure these aromatic 
sequences adopt as part of the viral fusion proteins. Impor- 
tantly, the pattern of hydrophobic/hydrophilic amino acids 
in the membrane proximal aromatic domains of class I viral 
fusion proteins is not consistent with the formation of 
amphipathic a helices, a motif commonly found in membrane- 
interacting helices. Instead, the aromatic residues in these 
aromatic domains appear to have a preference for an 
alternating or dyad repeat pattern (53), which is consistent 
with a extended conformation in the interfacially bound state. 
Our data for the CoV aromatic peptides are consistent with 
the idea that these aromatic sequences can have biologically 
relevant membrane activity without a significant degree of 
regular secondary structure. 

A working model for the process of SARS CoV virion/ 
cell fusion can be extrapolated from this work and from work 
in several other RNA virus systems (J0, 22, 24, 29, 42, 46, 
$2). After the fusion protein binds to the cellular receptor, 
the fusion protein undergoes a series of rearrangements, 
including the release of the receptor-binding subunit (S1) 
and trimerization and extension of the N helices of the S2 
subunit toward the target membrane. The next critical step 
in the fusion process is the insertion of the fusion peptide 
into the target cell membrane (83). The fusion peptide of 
any CoV has not been conclusively identified. Weiss and 
co-workers identified three sequences of MHV with clusters 
of hydrophobic amino acids that they postulated may 
represent the fusion peptide(s) of this model CoV (26). The 
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FiGuRE 7: Hypothetical mechanism for the involvement of the 
aromatic domain in SARS CoV virion/cell fusion. Class I viral 
fusion proteins have a fusion peptide (red) at the amino terminus 
and two extended © helices (N helix, orange; C helix, yellow) and 
most have an aromatic-rich domain (green) proximal to the TM 
anchor (indigo). After receptor binding with the S1 subunit, the 
helical domains of S2 “snap back”, bringing the viral and cell 
membrane in closer proximity. The rearrangement of the $2 protein 
into the six-helix bundle confirmation results in the alignment of 
the fusion peptide, aromatic domain, and TM anchor, constituting 
a contiguous track of hydrophobic sequences that can facilitate the 
flow of lipid between the two membranes and the eventual fusion 
of the viral and host cell membranes. 


WW interfacial hydrophobicity plot of the SARS CoV S82 
subunit (Figure 2A) shows a region of 19 amino acids (amino 
acids 770—788) with a score of 3.07 kcal/mol. On the basis 
of its similarities to the known HIV fusion peptide, we 
hypothesize that this region could potentially function as the 
SARS CoV fusion peptide. 


After insertion of the fusion peptide into the target cell 
membrane there is a rearrangement of the fusion protein, 
whereby a trimer-of-hairpins is formed (2/, 23). Subse- 
quently, the apposing membranes are brought close enough 
together by the assembly of a six-helix bundle and spontane- 
ously undergo first hemifusion and then complete fusion. 
The resulting structure is believed to align the fusion peptide, 
aromatic domain, and TM anchor (30), resulting in viral/ 
cell membrane fusion. On the basis of the observations that 
(1) the aromatic peptides of coronaviruses and other viruses 
(10, 55, 56, 58) strongly interact with and disrupt membranes 
and (2) the overall conservation of these sequences in many 
viruses, we hypothesize that the aromatic domain of class I 
viral fusion proteins plays a critical role in viral/cell 
membrane fusion. Specifically, we propose that the aromatic 
domain, in conjunction with the fusion peptide and the TM 
anchor, forms a continuous track of hydrophobic, membrane- 
interacting surfaces that provide a low-energy (low-barrier) 
path for lipid flow and membrane fusion during virion/cell 
fusion (Figure 7). 
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